**Volume Title** 

ASP Conference Series, Vol. **Volume Number** 
**Author** 

© **Copyright Year** Astronomical Society of the Pacific 



Predicting the Evolutionary Descendents of Sequence E Stars 

J. D. Nie 1 - 2 , P. R. Wood 1 , and C. P. Nicholls 1 

1 Research School of Astronomy and Astrophysics, Australian National 
University, Cotter Road, Weston Creek ACT 2611, Australia 

^Department of Astronomy, Beijing Normal University, Beijing, 100875, China 

Abstract. Sequence E variables are close binary red giants that show ellipsoidal light 
variations. They are likely to terminate their red giant evolution by a common envelope 
(CE) event when the red giant fills its Roche lobe, and produce close binary Planetary 
Nebulae (PNe). We made a Monte Carlo simulation to predict the fraction of Planetary 
Nebulae Nuclei (PNNe) that are post-CE binaries, using the observed frequency of 
sequence E binaries in the LMC to normalize our calculations. We find that 10-16% of 
PNNe should be short period, post-CE binaries. 



1. Introduction 



Sequenc e E variables are ellipsoidal binaries wh ich follow a loose period-luminosity 
relation dWood et al.lll999l : ISoszvnski et al.ll2004l) . In the LMC, they make up 0.5 to 2% 
of luminous red giant stars. Most of the sequen ce E stars lie on the first giant branch 
in the case of low mass stars ( Wood et apl999h . and some of them also evolve to the 
AGB dSoszynski et al.l 120041 : ISoszynskl 20071) . While all the sequence E stars show 
ellipsoidal varia tions, about 7% of t hem show eclipses in addition to the ellipsoidal 
light variations dSoszynski et al.ll2004 ). In the ellipsoidal binary systems, the primary 
star is a red giant substantially filling its Roche lobe, and the secondary star is an unseen 
companion orbiting the primary star. Due to the tidal interaction, the primary star is 
distorted by the companion, causing an ellipsoid or pear-like shape. As the system 
orbits around, the change in the apparent surface area gives rise to the ellipsoidal light 
variation. 

Sequence E stars are likely the precursors of close binary Planetary Nebulae (PNe). 
In the ellipsoidal binary systems, the red giant primary star has already substantially 
filled its Roche lobe. As the star evolves, it will experience a Roche lobe overflow and 
then faces a common envelope (CE) event. The primary star will terminate its red giant 
phase by ejecting the CE and it will produce a close binary PN. 

At the present tim e the most favoured theory explaining the nonspherical PNe is 
the bi nary hypothesis dBond & Livid 1 1 990l: lYungelson et all 1 19931 : ISokerlll997l : iBondl 
20001 ; IZiilstral 120071 : Ide Marcdl2009h . Theoietical considerations suggest that a binary 
with CE will eject the entire envelope to produce an asymmetrical PN, with elliptical or 
bipol ar shape. In fact , PNe with close binary nuclei have been confirmed by ob serva- 
tions dBond & Grauer!ll987l : lBond et aLlll992i^BonJl994l2000l : lMiszalski et al.ll2009h . 
However, about their fraction, little is known. Bond and collaborators found about 10- 
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15% of PNe with close binary central stars, while iMiszalski et "ail (2009) found 12-21% 
of Planetary Nebulae Nuclei (PNNe) are close post-CE binaries. 

Sequence E stars are likely to prod uce close binary PNe, and their fraction is 
known. In the LMC, IWood et all (l999) found 0.5% of the red gian ts on the top 1 



magnitude of the RGB in the LMC are sequence E variables. Similarly. iSoszyhski et al. 



(2004) and lSoszvnskl d2007l) found 1-2% (we use 1.5%) of red giants in the LMC show 
ellipsoidal variations. Moreover, for these fractions, we found the detectability limit of 
the full light curve amplitude is ~0.05 mag for Macho Red band (Mr), and ~0.025 mag 
for OGLE / band. It is our aim to use the observed fraction of sequence E stars among 
all the red giants to estimate the fraction of PNe with close binary central stars that will 
be produced by a CE event. 



2. Monte Carlo Simulation 

A Monte Carlo simulation is made to predict the fraction of close binary PNe. One 
million red giant binaries are initially generated by using the observed orbital elements 
distributions. All the binaries are evolved up to the RGB and AGB and their evolu- 
tionary fates are examined to find which stars would produce close binary PNe. The 
fraction of close binary PNe is estimated by using the observed fraction of sequence E 
stars on the top 1 magnitude of the RGB. 



2.1. Generating the binary systems 



The simulation requires as input the orbital elements distributions of the binary sys- 
tems. However, these are poorly known in the LMC, so we adopt the distributions from 
binaries in the solar vicinity. Input from LMC sources are used when available. 

(1) We consider the Initial Mass Function (IMF) as well as the star formation 
history to get the initi al mass distribution of the primary star. The IMF is assumed to 
follow ISalpeterl (Il955l) 's power law. The LMC star formation history is adopted from 
iBertelli et al. (1992): a star burst begins ~4 Gyrs ago and ceases ~0.5 Gyrs ago, with the 
ratio of burst to quiescent st ar formation rate being 10. According to the evolutionary 
tracks of iGirardi et al.1 ( 20001) . a 4 Gyrs old star has a mass of 1.3 M©, and a 0.5 Gyrs 
old star has a mass of 3.0 M G . The mass range is set from 0.9 M Q to 1.85 M Q for RGB 
stars, where 0.9 M© is th e mass for red giants born in the early universe of age 13.7 
Gyrs (ISpergel et al. 12003 ). and 1.85 M Q is the upper limiting mass for red giants with 
electron degenerate helium cores on the first giant branch. For AGB stars, a full mass 
range from 0.9 to 3.0 M Q is allowed. 



(2) The mass ratio of the binary systems is drawn from the distribution of lDuquennoy & Mayor 

(Il99ll) . which follows a Gaussian-type relation, with its peak at 0.23. 



(3) The orbital period of the binary systems is drawn from the distribution of 

which follows a Gaussian-type relation, with its peak at 



Duquennoy & Mayo: 



period o 



173 years. 

(4) We assume the eccentricity is zero. 

(5) The orbital separation is calculated by using the equation of the orbital motion. 

(6) The orbital inclination is obtained assuming a random orientation of the orbital 

pole. 
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2.2. The properties of the binary systems 
2.2.1. Mass loss 

We consider mass loss from the RGB stars via a stellar wind. We use the empirical for- 



mulation by lReimersI (11975!) to calculate the mass loss rate, but with the rate multiplied 



by a p arameter n which is set equal to 0.33 (llben & Renzinilll983l : iLebzelter & Wood 



2005). In order to know the amount of mass lost per magnitude of evolution up the red 
giant bran ch, we need to ca lculate the evolution rate. According to the evolutionary 
track of Girar di et ail ((2000), the evolution rate for RGB stars is set to dM bo i/d? » 0.15 



mag/Myr. Due to the mass loss from the primary red giant, the binary system loses 
its orbital angul ar momentum. The orbital evolution of the system is calculated with 
equation (20) in lHurlev et"aL1 d2002h . 



AGB mass loss driven by the superwind is not considered explicitly. At the AGB 
tip, stars will lose their mass by a superwind, causing the termination of the AGB 
evolution by ejecting the envelope, just as happens for single stars, producing single 
PNe or wide binary PNe. We treat the AGB superwind simply as terminating AGB 
evolution at a specific luminosity. 

2.2.2. Stellar radius 

For an ellipsoidal binary system containing a red giant, the requirement for the detection 
of the ellipsoidal variability is the minimum fractional filling of the Roche lobe. In 
order to know the minimum filling factor of the Roche lobe, we use the light curve 
generator nightfaliQ to model the observed ellipsoidal light variations of partial-Roche 
lobe filling systems as a function of binary parameters. The input binary parameters 
are typical of sequence E stars. We find the relation between the light curve amplitude 
(Mr) and Roche lobe filling factor '/' i s: 

AM R = (0.221/ 4 + 0.005) x (1.44956<? - 25 - 0.44956) x sin 2 i , (1) 

for 0.5 < f < 0.9, 0.1 < q < 1.5, < i < |. We also find the relation between AM R 
and AI is well approximated by: 

A7 = 0.87xAM s , (2) 

where AI is the full light curve amplitude in the I band. 

Given the minimum amplitude for detectable light variations (0.05 mag in the Mr 
band or 0.025 in the I band), the minimum radius filling factor / m ; n of a red giant 
in a binary system with detectable ellipsoidal light variation can be determined from 
equations (O and ©. With minimum Roche lobe filling factor / m j n , the just-detectable 
stellar radius is express as: 

-^min — X f m { n , (3) 



where 7?l is the equivalent radius of the Roche lobe (Eggleton 1983). When the Roche 
lobe is filled, the maximum stellar radius is: 

^max = ■ (4) 



'http://www.hs.uni-hamburg.de/DE/Ins/Per/Wichmann/Nightfall.html 
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2.2.3. Effective temperature and luminosity 

The effective temperature and luminosity of the primary red giants are calculated sep- 
arately for O-rich and C-rich stars. To obtain these, we use the intermediate age LMC 
globular cluster NGC 1978 as a template for LMC red giant properties. 



For O-rich stars, the giant branch s lope is obtained by u sing the data of lKamath et al 



( 20101) . With the evolutionary tracks of lGirardi" et all (12000). we make a mass correction 



for the effective temperature. The zero point of the effective temperature r e ff (Mt, i, m) 
is estimated by fitting the observed HR diagra m of sequence E stars . For C-rich stars, 
we use the giant branch slope obtained from iKamath et all d2010l) . The zero point 
of the r e ff (Mboi, m) relation for C-rich stars is obtained by equating the temperatures 
of the O-rich and C-rich stars at the transition luminosity from t he O-rich to C-rich 
star s. The transition luminosity is calculated by using the data of iFrogel et all (1 19901) 



and IVassiliadis & Wood! dl994h . with the distant modulus of the LMC equal to 18.54 
dKeller & Woodl2006h . 



With L = Ano-R 2 T A eS , logL/L Q = -0.4(M bo i - 4.75) and the T eS (M hoh m) from 
above, we derive the bolometric magnitude for O-rich and C-rich stars as a function of 
m and R. If we substitute R with R m i n or R max , then we get the minimum or maximum 
bolometric magnitude of the red giants with just-detectable light variations or with 
Roche lobes that are just full, respectively. 



2.3. Scenario for finding close binary PNe 

An ellipsoidal variable almost filling its Roche lobe will experience a Roche lobe over- 
flow as the star evolves. As a result of the Roche lobe overflow, the binary system will 
suffer a CE event, leading to the ejection of the entire envelope. A close binary PN will 
be produced if the primary star is on the AGB, or a post-RGB binary system will be 
formed if the primary star is on the RGB0. For RGB binaries, since the evolution is too 
slow after envelope ejection, the system will be unable evolve to high effective temper- 
ature quickly enough to ionize the ejected envelope before it disperses. On the other 
hand, post-RGB stars will evolve slowly through the instability strip where they will be 
seen as Population II Cepheids or RV Tauris. We note that some Populati on II Cepheid s 
and RV Tauris do have luminosities below the RGB tip luminosity (lAlcock et al.ll 1998 ). 



3. Results 

With the model described in Sect. |2J we calculated the fraction of PNNe that are close 
binaries. The fraction of sequence E stars in the LMC was used to normalize the calcu- 
lation. In order to get the observed fraction of sequence E stars, we added single stars 
until the measured fraction is equal to the observed. 

The predicted fraction of close binary PNe is -10 or 16%, depending on whether 
we use Wood's (0.5%) or Soszyhski's (1.5%) frequency. It indicates that short-period 
post-CE binaries are a small fraction of PNe central stars. 

To test the simulation, we predicted the distributions of the period and velocity 
amplitude of sequence E stars and compared them with the observations. Fig. [U(left 
panel) is the period distribution from the model and observation (OGLE II). It shows 



2 We did not consider the coalescence of the two companions. 
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Figure 1 . The distributions of the orbital period (left) and full velocity amplitude 
(right) of sequen ce E stars. Solid line d enotes distribution from the observation (left: 
OGLE II; right: Nicho llset al.1 ((2010)), dotted line denotes distribution from the 
model by using Wood's frequency, short dashed line denotes distribution from the 
model by using Soszyhski's frequency. 



that the simulation reproduces the observation well. The period ranges from ~30 to 
1000 days, with the peak located at -250 days. We also predicted the full velocity 
amplitude distribution (right panel in Fig. [[J. The distribution from the model shows 
a good consistency with the observational results of iNicholls et all (|2010|) although the 
observed numbers are small. 

An updated and improved version of the results presented here is currently being 
prepared for publication (Nie et al. 201 1, in preparation). 
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